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ABSTRACT 

We present a new method how to discriminate the matter content of parsec-scale jets 
of active galactic nuclei. By constraining the kinetic luminosity of a jet from the observed 
core size at a single very long baseline interferometry frequency, we can infer the electron 
density of a radio-emitting component as a function of the composition. Comparing 
this density with that obtained from the theory of synchrotron self-absorption, we can 
determine the composition. We apply this procedure to the five components in the 3C 345 
jet and find that they are likely pair-plasma dominated at 11 epochs out of the total 21 
epochs, provided that the bulk Lorentz factor is less than 15 throughout the jet. We also 
investigate the composition of the 3C 279 jet and demonstrate that its two components 
are likely pair-plasma dominated at three epochs out of four epochs, provided that their 
Doppler factors are less than 10, which are consistent with observations. The conclusions 
do not depend on the lower cutoff energy of radiating particles. 

Subject headings: galaxies: active — quasars: individual (3C 345) — quasars: individual 
(3C 279) — radio continuum: galaxies 
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1. Introduction 

The study of extragalactic jets on parsec scales is 
astrophysically interesting in the context of the ac- 
tivities of the central engines of active galactic nuclei 
(AGN). In particular, a determination of their mat- 
ter content would be an important step in the study 
of jet formation, propagation and emission. There 
are two main candidates for their matter content: 
A 'normal plasma' consisting of (relativistic or non- 
relativistic) protons and relativistic electrons (Celotti 
and Fabian 1993; see also Gomez et al. 1993, 1994a, b 
for numerical simulations of shock fronts in such a 
jet) and a 'pair plasma' consisting only of relativis- 
tic electrons and positrons (e.g., Kino and Takahara 
2004). Discriminating between these possibilities is 
crucial for understanding the physical processes oc- 
curring close to the central engine (presumably a su- 
pcrmassive black hole) in the nucleus. 

Very long baseline interfcrometry (VLBI) is uniquely 
suited to the study of the matter content of pc-scale 
jets, because other observational techniques cannot 
image at milliarcsecond resolution and must resort to 
indirect means of studying the active nucleus. So far, 
there have been two observational approaches with 
VLBI to discriminate jet compositions: polarization 
(Wardle et al. 1998) and spectroscopic (Reynolds et 
al. 1996) techniques. 

Wardle et al. (1998) observed the 3C 279 jet at 15 
and 22 GHz. They examined the circular polarization 
(<~ 1%) due to Faraday conversion of linear to circu- 
lar polarization caused by low-energy electrons and 
derived that the lowest Lorentz factor, 7 m i n , of radi- 
ating particles, should be much less than 10 2 . They 
also analyzed the linear polarization (<~ 10%), which 
strongly limits internal Faraday rotation to derive 
that 7 m ; n ^> 10 2 must be held for a normal-plasma 
composition and that no constraint would be obtained 
for a pair-plasma composition. (Note that an equal 
mixture of electrons and positrons can produce Fara- 
day conversion, but not rotation.) From these argu- 
ments, they concluded that component CW of the 
3C 279 jet consists of a pair plasma with 7 m j n <C 10 2 . 

On the other hand, Reynolds et al. (1996) ana- 
lyzed historical VLBI data of the M87 jet at 5 GHz 
(Pauliny-Toth et al. 1981) and concluded that the 
core is probably dominated by an e* 1 plasma. In 
the analysis, they utilized the standard theory of syn- 
chrotron self absorption (SSA) to constrain the mag- 
netic field, B [G], and gave B = 0.2 G (line B in fig- 



ure 1). They also considered the condition that the re- 
solved core becomes optically thick for self-absorption 
and derived another constraint on TV* and B. Using 
the observed angular diameter of the core (0.7mas), 
the total core flux density (1.0 Jy) at 5 GHz, and the 
reported viewing angle (30° — 40°), and assuming that 
the spectral index is 0.5 (see end of this section), they 
derived N*B 2 > 2<5 rn 1 { in cgs unit (line C in figure 1), 
where <5 m ax (= 2 for the M87 core) refers to the upper 
limit of the Doppler factor of the fluid's bulk motion. 

This condition is, however, applicable only for the 
VLBI observations of M87 core at epochs September 
1972 and March 1973. Therefore, in order to apply 
the analogous method to other AGN jets or to the 
M87 jet at other epochs, we must derive a more gen- 
eral condition. 

On these grounds, Hirotani et al. (1999, hereafter 
Paper I) generalized the condition N* B 2 > 2S~ 2 K and 
applied it to the 3C 279 jet on parsec scales. In that 
paper, they revealed that the core and components 
C3 and C4, of which spectra are reported, are likely 
dominated by a pair plasma. It is interesting to note 
that the same conclusion was derived by an indepen- 
dent method by Wardle et al. (1998) on component 
CW in the 3C 279. 

Subsequently, Hirotani et al. (2000, hereafter Pa- 
per II) applied the same method to the 3C 345 
jet. Deducing the kinetic luminosity, Lkin, from a 
reported core-position offset (Lobanov 1998), they 
demonstrated that components C2, C3, and C4 at 
epoch 1982.0, C5 at 1990.55, and C7 at four epochs 
are likely dominated by a pair plasma. In the present 
paper, we propose a new scheme to deduce the kinetic 
luminosity of an unresolved core and re-examine the 
composition of the two blazers 3C 345 and 3C 279, 
adding spectral information at 11 epochs for the 
3C 345 jet to Paper II and and 2 epochs for the 3C 279 
jet to Paper I. 

In the next section, we give a general expression of 
lines B and C in figure 1. We then describe in § 3 a 
new method how to infer Lkin from the core size ob- 
served at a single VLBI frequency; the inferred Lkin 
gives lines Di and D2 in figure 1, depending on the 
composition. Once Lkin is obtained for the core by 
this method, we can apply it to individual jet com- 
ponents, assuming a constant ikin along the jet. In 
sections 4 and 5, we examine the two blazers 3C 345 
and 3C 279. In the final section, we discuss the pos- 
sibility of the entrainment of the ambient matter as a 
jet propagates downstream. 
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We use a Hubble constant Hq = 65ft. km/s/Mpc 
and go = 0.5 throughout this paper. Spectral index 
a is defined such that S v oc v +a . 

2. Synchrotron self-absorption constraint in 
the jet region 

In this paper, we model a jet component as a homo- 
geneous sphere of angular diameter 9d, containing a 
tangled magnetic field B [G] and relativistic electrons, 
which give a synchrotron spectrum with optically thin 
index a. In § 2.1, we constrain B (i.e., line B in fig. 1) 
for arbitrary optical thickness by the theory of syn- 
chrotron self-absorption. Once B is obtained, we can 
compute the electron number density A^*(SSA) (i.e., 
line C) as will be described in § 2.2 with the aid of the 
computed optical depth at the turnover frequency. 

2.1. Magnetic field strength 

The magnetic field strength, B, can be computed 
from the synchrotron spectrum, which is character- 
ized by the peak frequency, f m , the peak flux density, 
5 m , and the spectral index, a, in the optically thin 
regime. Marscher (1983) considered a uniform spher- 
ical synchrotron source and related B with z/ m , S m , 
and a for optically thin cases (0 > a > —1.0). Later, 
Cohen (1985) considered a uniform slab of plasma as 
the synchrotron source and derived smaller values of 
B for the same set of (v,S m ,a). In this section, we ex- 
press B (i.e., line B in fig. 1) in terms of (v,S m ,a) for 
a uniform sphere for arbitrary optical thickness (i.e., 
for arbitrary a < 0). We assume that the magnetic 
field is uniform and that the distribution of radiating 
particles are uniform in the configuration space and 
are isotropic and power-law (eq.[19]) in the momen- 
tum space. 

For a uniform B and TV*, the transfer equation 
gives the specific intensity 



L * =A „*5/2 [l_ cxp (_ ai/a;o *)] 



(1) 



where 
A(a) 



3\ °e a(a) 
2 



c C(a) \27TTO e c 



-3/2 



B- 1 / 2 , (2) 



and xo* gives the physical thickness of the emitting 
region along the line of sight; e, m c , and c refer to the 
charge on an electron, the rest mass of an electron, 
and the speed of light, respectively. The coefficients 
a(a) and C(a) are given in table 1; a quantity with 



an asterisk is measured in the co-moving frame, while 
that without an asterisk in the observer's frame. At 
the distance 9 away from the cloud center (fig. 2), 
the fractional thickness, which is Lorentz invariant, 
becomes 



||- = C os(0 + O = 



sin 9 



sin(0 d /2) 



(3) 



where R* cos£+ (R* sin£/ svaO) cos 9 = R* / sin(6 l d/2) 
is used in the second equality; 6>d is the angular diam- 
eter of the component in the perpendicular direction 
of the jet propagation. Thus, the specific intensity at 
angle 9 away from the center, becomes 



1 + z 



1/2 



Av 5 ' 2 



1 — exp 



-r„(0)Wl- 



sin( 



sin(0 d /2) 



(4) 



where r„(0) = a u * ■ 2R* is the optical depth for 
9 = (or b = 0). Averaging over pitch angles of the 
isotropic electron power-law distribution (eq. [19]), 
we can write down the absorption coefficient in the 
co-moving frame as (Le Roux 1961, Ginzburg & Sy- 
rovatskii 1965) 



C(a)r 2 k, 



v 



v* 

2//Vv, „2 



(-2a+3)/2 



(5) 



where v Q = c/r Q = c/[e 2 /(m c c 2 )] = 1.063 x 10 14 GHz 
and v-q = eB/(2irm e c). In deriving equation (4), we 
utilized Ivjv 3 = I v * jv* 3 and vjv* = 6/(l + z), where 
z represents the redshift and the Doppler factor, <5, is 
defined by 

6 = r(l-/3cos^)' (6) 



where T = 1/ y/l — {3 2 is the bulk Lorentz factor of 
the jet component moving with velocity /3c, and ip 
refers to the viewing angle. 

Even if special relativistic effects are important, 
the observer find the shape to be circular. We thus 
integrate I u (9) cos 9 over the emitting solid angles 
2n sin 9d9 in < 9 < 9^/2 to obtain the flux den- 
sity 



S v 



7rsm 



Sd/2 



I„{6) cos 9 sin 9d9 
1 + z 



3 



-Tv(0)VI=X 



dx, (?) 



where x = [sin 9/ sm(6> d /2)] 2 . 

Differentiating equation (7) with respect to v, and 
putting dS u /dv to be 0, we obtain the equation that 
relates r„(0) and a at the turnover frequency, z/ m , 



-t„(0)v/I=x 



dm 



We denote the solution t u (0) at i/ = v m as r m (0), 
which is presented as a function of a in table 1. 

It is worth comparing equation (7) with that for a 
uniform slab of plasma with physical thickness R* . If 
the slab extends over a solid angle 7r(0.5^d/rad) 2 , the 
flux density becomes 



S,, 



4 V rad 



1 + z 



1/2 



(9) 

The optical depth t„(0) = a* u -2R* does not depend on 
for the slab geometry. Therefore, comparing equa- 
tions (7) and (9) we can define the effective optical 
depth, (t„), for a uniform sphere of plasma as 



-<r„> = 



dx- 



(10) 



We denote (tv) at v — z/ m as (r m ), which is presented 
in table 1. Note that (r m ) is less than T m (0) due to 
the geometrical factor y/l — %. 

Using r m (0), we can evaluate the peak flux density, 
S m , at ^ = z/ m and inversely solve equation (7) for B. 
We thus obtain 



b = 10-^) (-*lv (— v a — 



where #d -C 1 rad is used and 



z 
(11) 



6(a) = 3.98 x 10 3 



-la r 



a(a) 



n 2 



C(«) 



1 . 



1 - e 



d X . (12) 



The values are tabulated in table 1; they are, in fact, 
close to the values obtained by Cohen (1985), who 
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Fig. 1. — Constraints on (B,N*) plane imposed by 
synchrotron self-absorption and total kinetic luminos- 
ity considerations (see fig. 1 in Reynolds et al. 1996). 
The VLBI surface brightness constraint on B gives 
line B, while the condition that the a jet component 
(or the resolved core in Reynolds et al. 1996) be- 
comes optically thick to self- absorption gives line C. 
The horizontal lines are the proper densities derived 
from the kinetic luminosity; line Di and D2 corre- 
spond to a pure pair and a pure normal plasma, re- 
spectively. 




Fig. 2.- 
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Flux emitted from a spherical cloud of ra- 
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presented b(a) for a slab geometry with —a < 1.25. 
This is because the averaged optical depth, (r„) at 
v = Vm, for a spherical geometry becomes compara- 
ble with that for a slab geometry (Scott & Readhead 
1977). 

We could expand the integrant in equation (7) in 
the optically thin limit r m (0) <C 1 as 



1 



!(0)VT 



x, 



(13) 



to obtain 6(-0.5) = 3.34 and fc(-l.O) = 3.85, for ex- 
ample. This optically thin limit (r m <C 1) was consid- 
ered by Marscher (1983), who gave 6(a) for a > —1.0 
(i.e., for r m (0) < 1.0). 

2.2. Electron density 

We next consider how to constrain N* (i.e., line C 
in fig. 1) from the theory of synchrotron self-absorption. 
To examine the lower limit of iV*(SSA), we consider 
a conical jet geometry (fig. 3) in this section, even 
though we apply the results (eqs. [22]- [24]) to indi- 
vidual jet components, which would be more or less 
spherical rather than conical. The optical depth r 
for synchrotron self absorption at distance p from the 
injection point, is given by 



psmx 



psmx 



i(tp + x) sin(ip - x) 



(14) 



where p> is the viewing angle and a v [1/cm] refers to 
the effective absorption (i.e., absorption minus stim- 
ulated emission) coefficient. For a small half opening 
angle (x <C 1), this equation can be approximated as 



T v (p) = 2— a,, 

sin ip 



(15) 



where x/sin^J denotes the projected, observed open- 
ing angle. Noting that psinx represents the trans- 
verse thickness of the jet, we find that equation (15) 
holds not only for a conical geometry but also for a 
cylindrical one. 

Since r„ and psinx are Lorentz invariants, we ob- 
tain 

= (i 6 ) 
sin ip sin p>* 

Since va v is also Lorentz invariant, equation (16) 
gives 

sin p* v 5 
sin p v* 1 + z 



(17) 



Combining equations (15) and (17), we obtain 

1 + z 2p sin x * _ 1 + z 8 d D L t 
5 sin^ v 5 s'mp(l + z) 2 VJ 

(18) 

where D^/(l + z) 2 is the angular diameter distance 
to the AGN; 9^ is measured in radian unit (i.e., not 
in milliarcsecond) . Equation (18) holds for a coni- 
cal geometry with an arbitrary small opening angle, 
including a cylindrical case. 

We assume that the electron energy distribution is 
represented by a power-law, 



dm 



= Kr 



2a- 1 



(7min < 7 < 7 max 



where N* refers to the proper electron number den- 
sity. Integrating dN^/dj from 7 m i n to j max , and as- 
suming 7 m ax 3> 7min and a < 0, we obtain 



/V* — ' mm ;„* 



(20) 



The coefficient C(a) is given in Tabic 1 of Gould 
(1979) and also in table 1 in this paper. 

Substituting equation (5) into (18), and assuming 
7min < 7max, we obtain 



N * B -<*+l.5 



m e c I e 
e 2 \ 2ttiti c 



-1.5+a 



2a 



C(a) -2a 



sinp (1 + z) 2 f 1 + z 



-a+2.5 



(21) 



Evaluating v at v m , and combining with equation 
(11), we obtain AT*(SSA) in cgs unit as (see also 
Marscher 1983 for relatively optically thin cases, —a < 



injection 
piont 




Fig. 3. — Schematic figure of a conical jet in the core 
region with half opening angle x m t ne observer's 
frame. 



5 



Table 1: Tabic of constants 



a 





-0.25 


-0.50 


-0.75 


-1.00 


-1.25 


-1.50 


-1.75 


-2.00 


a 


0.2833 


0.149 


0.103 


0.0831 


0.0740 


0.0711 


0.0725 


0.0776 


0.0865 


C 


1.191 


1.23 


1.39 


1.67 


2.09 


2.72 


3.67 


5.09 


7.23 


7- m (0) 





0.252 


0.480 


0.687 


0.878 


1.055 


1.220 


1.374 


1.519 







0.167 


0.314 


0.445 


0.564 


0.672 


0.770 


0.862 


0.946 


b 





1.66 


2.36 


2.34 


2.08 


1.78 


1.51 


1.27 


1.08 



1.0) 



iV*(SSA) = e(a) 

x vghz/ VJy 



7 min sm ip / & d 



—2a c?a \mas / 

2a+3 / r \ 2a-3 



1 + Z 



where 



d A 



zqo + {qa - 1)(-1 + V'Mqz + 1) 



h(l + z) 2 q 2 



(22) 



(23) 



e(a) = 1.71 x 10~ 9 x [2.79 x 10" 8 6(a)] 



a-1.5 



C(a) 
(24) 

If the jet component is discontinuous (e.g., spherical), 
equation (22) gives the lower limit of -ZV*(SSA), be- 
cause the path length becomes smaller than 2p sin \l sin ip 
(eq. [15]). 

It looks like from equation (22) that 7V*(SSA) de- 
pends on #d and v m very strongly. In the case of 
a = —0.75, for instance, we obtain 7V*(SSA) oc 



d 



~io, 



Nevertheless, as a radio-emitting com- 



ponents evolves along the jet, its #d increases due 
to expansion while its v m decreases due to (syn- 
chrotron+adiabatic) cooling. As a result, these two 
effects partially cancel each other and suppress the 
variations of A"*(SSA) along the jet. 

3. Kinetic luminosity in the core region 

In this section, we present a new method how to 
deduce the kinetic luminosity, Lkin, of a jet in the un- 
resolved VLBI core from the perpendicular core size, 
utilizing the core-position offset among different fre- 
quencies as an intermediate variable in manipulation. 
Once Lkin is obtained for the core, we can apply it to 
derive lines Di and D2 for individual jet components, 
assuming a stationary jet ejection with constant Lkin- 

3.1. Scaling Law 

Let us introduce a dimensionless variable r = p/r\ 
to measure p in parsec units, where n = 1 pc. We 



assume that the electron number density scales with 
r as 



A c * = N\r~ 



(25) 



where N\ refers to the value of TV* at r — 1 (i.e., at 1 
pc from the core in the observer's frame). If particles 
are neither created nor annihilated in a conical jet, 
conservation of particles gives n = 2. In the same 
manner, we assume the following scaling law: 



B = B-,r~ 



(26) 



where B\ refers to the values of B at r = 1. In a 
stationary super-fast-magnetosonic jet, the magnetic 
flux conservation gives m = 1. Note that we assume 
such scaling laws only in the unresolved VLBI core, 
not in the jet region, for which we consider iV*(SSA). 

We further introduce the following dimensionless 
variables: 



2 A 1; 
1 

2-Km e c vq 



Xn = ri r. 
eB 1 



(27) 



Utilizing equation (5), and rewriting k* and isq in 
terms of in and xb, we obtain from the left equality 
in equation (18) 



C(a)- 



2X 



-2a 



Sm (fi 7 m in 2Q 



1 



(28) 



where e = 3/2 — a. Note that both and xq should 
be distinguished from N*(SSA) and B presented in 
section 2, because we are considering the core region 
in this section. 

At a given frequency v 1 the flux density will peak 
at the position where t v becomes unity. Thus setting 
t = 1 and solving equation (28) for r, we obtain the 
distance from the VLBI core observed at frequency v 
from the injection point as (Lobanov 1998) 

l/fer 



r{v)= x B kb f(a)-^ 



(29) 
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where 
/(«) 



C{a) 



2x 



-2a 



sm ip 7 min 
3 



2a 



1 + Z 



I N 



2a 



5 — 2a 
(3 - 2a)m + 2n - 2 



3.2. 



5- 2a 
Core-Position Offset 



(30) 
(31) 

(32) 



If we measure r(v) at two different frequencies (say 
f a and Vb), equation (29) gives the dimensionless, pro- 
jected distance of r(z/ a ) — r{vb) as 



Ar D 



= (V b ./Vo) 



r(f b )] sint^ 

1/fcr 



\Vfcr 



l/fer 



l/fe r l/fe r 

fa I/ b 



■sin (^.(33) 



Here, Ar pro j is in pc units; therefore, riAr pro j repre- 
sents the projected distance of the two VLBI cores (in 
cm, say). That is, r\ Ar pro j equals 4.85T0 _9 Ar mas DL/ (1+ 
z) 2 in equation (4) in Lobanov (1998), where Ar mas 
refers to the core-position difference in mas. Defining 
the core-position offset as 



Q rv = nAr, 



l/fer 1/fcr 
V b 

P r °J l/fe r 1/fcr 



we obtain 



a 

n 



^ = (4v^) 1 / fer si 



sm ip 



(34) 



(35) 



Setting Vb — ► oo in equation (33), we can express 
the absolute distance of the VLBI core measured at 
v from the central engine as (Lobanov 1998) 



= (") = 



-l/fer 



(36) 



n sm ip 

To express ib in terms of in and Q rv , we solve equa- 
tion (35) for ib to obtain 

*b= -^-) (.f-o)- 1/feb . (37) 

Note that in is included in / = f(a). 

We next represent in and ib (or equivalcntly, N\ 
and Bi) as a function of Q rv . To this end, we pa- 
rameterize the energy ratio between particles and the 
magnetic field as 



B 

N*^ min m c c 2 = K—. 

an 



(38) 



When an energy equipartition between the radiating 
particles and the magnetic field holds, we obtain 



K = 



7n 



(7-> 



(39) 



where the averaged electron Lorentz factor, (7_) be- 
comes 



(7-) = 



f 

Jin 



7 • fc e *7 2a -^7 



7mi 

2a 



2a + 1 



fc e *7 2a - 1 d 7 

(7max/7min) 2a "^ 1 — 1 
V /max/ /mm ) 1 



7n 



.(40) 



for a < 0. We assume a rough energy equipartition 
holds and adopt 



K 



2a + 1 ( 

Trnax/ Tmin) 



l 



2a f-v /'v • ^a+i _ i ' 
^a ^ [max/ /miny 1 



(41) 



In the limit a 



-0.5, the right-hand side tends to 



l/ln(7 max /7 min ), which is 0.1 if 7 n 



10 4.d4 7 



for instance. Note that we assume a rough energy 
equipartition in the core and not in the jet. 

In this paper, we assume that K/j m i D is constant 
for r. Then, equation (38) requires n — 2m. Sub- 
stituting N* = Nir~ 2m and B = Bir~ m into equa- 
tion (38), and replacing N\ and B\ with in and ib, 
we obtain 

7r K n 2 



I N = 



-IB 



2 7min 

Combining equations (37) and (42), we obtain 

-i (5-2a)/(7-2a) 



(42) 



I B 



irC(c 



n sm ip , 

x K n 



-2a 



sm ip 7 min r G 7 n 



l + z 



-2/(7-2a) 

(43) 



For ordinary values of a(< 0), ib decreases with in- 
creasing K, as expected. The particle number density, 
in, can be readily computed from equation (42). 

3.3. Kinetic Luminosity 

When the jet has a cross section ttR 2 ^ at a certain 
position, Lkin and N* are related by 

ikin = 7ri? ± 2 /3c-r7V*-(r-l) (( 7 -}m c c 2 + {l+)m+c 2 ) , 

(44) 
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where (7+) refers to the averaged Lorentz factors of 
positively charged particles, and m + designates the 
mass of the positive charge. 

Assuming that the particle number is conserved, 
we find that nR± 2 f3c ■ FN* is constant along a sta- 
tionary jet. Thus, provided that T > 1 holds, R± 2 N* 
can be replaced as 

^ 2 ^e* = ^iXi) 2 ^ = ^Xi 2 ^N, (45) 

1 1 To 

where and \i refer to T and x at 1 pc. Substituting 
equation (45) into (44), we obtain 



Ckin^/jr^r-i) (*I 

To 6 V X 



x_ 

v \ n sin ip 
nC{a) K n 



2(5-2a)/(7-2a) 



-2a 



Sin <p 7min To 7min 



2 a 



1 + z 



3/2- 



-4/(7-2a) 

(46) 



we adopt Cki n = tt 2 (l-)fn c c 3 /7 m i n for a pair plasma, 
while Ckin = 71-2 (7+) TO p c 3 /(27min) for a normal plasma, 
where to p refers to the rest mass of a proton. Equa- 
tion (46) holds not only for a conical jet but also for 
a general jet with T ^> 1 provided that the parti- 
cle number is conserved and that n = 2m holds. It 
should be noted that 7 m i n takes different values be- 
tween a pair and a normal plasma. 

3.4. Jet Opening Angle 

Let us now consider the half opening angle of the 
jet, and further reduce the expression of L^ n . As 
noted in § 3.3, we obtain 



X = 



R± = 0.5fld,corc(^ PL 
P " " " " P (1 + ^) 2 



(47) 



for a conical geometry, where #d,corc(^) is the angular 
diameter of the VLBI core at a frequency v in the 
perpendicular direction to the jet-propagation direc- 
tion on the projected plane (i.e., a VLBI map). The 
luminosity distance is given by 



D h = 4.61xlO^-V 1 goZ + (g0 - + ^1±A . 

9o 2 

(48) 

On the other hand, equation (36) gives 



Qrv —l/kr 

p = r corc ri = v 1 ' 

sm ip 



(49) 



Substituting equations (49), (48) into (47), we obtain 
1 \Q rv 



va ri sm ip 



11.1ft" 



.i?o* + («o - 1)(-1 + V'2<loz + 1) 



9o 2 (l + ^) 2 



'd,corc 

mas 



u l/kr 

"0 



(50) 



Since to = 1 is consistent with the theory of mag- 
netohydrodynamics, we will adopt k r = 1 in the rest 
of this paper. Then, it follows that the factor con- 
taining x^rvj simp in equation (46) can be simply 
expressed in terms of 6*d,core(^)- That is, to evalu- 
ate Lkin, we do not have to stop on the way at the 
computation of Q, rll , which requires more than two 
VLBI frequencies. We only need the VLBI-core size 
at a single frequency, 9 A _ cmc (v), m addition to T, <p, 
a, 7mi n , z, K, and Cki n - If we have to examine x, we 
additionally need f2 r!/ . The information of the com- 
position is included in Ckin- In the present paper, we 
assume Tx 2 = TiXi 2 - 

If k r deviates from unity, we have to independently 
examine \ and Q rl/ . In another word, we need at 
least two VLBI frequencies (e.g., 2 and 8 GHz) to 
measure the core-position offset in order to constrain 
ikin, even if we know the value of k r by some means. 

To constrain L in equation (46), let us briefly con- 
sider how to obtain its lower bound. If the appar- 
ent velocity of a radio-emitting component, /3 app = 
f3smip/(l — (3 cos ip), reflects the fluid velocity, the 
lower bound can be obtained by L m ; n = ^//3 app 2 + 1- 
On the other hand, the upper bound of T is not eas- 
ily constrained; thus, we must consider it source by 
source, using observations at different photon energies 
(e.g., X-rays) if necessary. 

3.5. Electron density deduced from kinetic 
luminosity 

Once Lkin is known by the method described above 
or by some other methods, we can compute the elec- 
tron number density as a function of the composi- 
tion. We assume that the i k in deduced for the core 
can be applied for the individual jet components. 
Then, solving equation (44), which holds not only 
in the core but also in the jet, for N*, and utiliz- 
ing R± = (6d/2)Di J /(l + z) 2 , we obtain for a normal 
plasma 



AC(nml) 



5.93 x 10 _2 ft' 



L 



eg y 

L(L - 1) Vmas / 
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x qo 2 (l + z) 2 y 

_zq + (qo ~ 1)(-1 + V' 2 qoz + 1). 

x — — — - ; \ ; cm~ 3 , (51) 

( 7+ ) 1 + (7-)m e /m p 

where L 4 6 = £ k in/(10 46 ergs s -1 ). If iV*(nml) be- 
comes less than the electron density deduced indepen- 
dently from the theory of synchrotron self-absorption 
(SSA), the possibility of a normal plasma dominance 
can be ruled out. 

It is worth comparing A*(nml) with TV* (pair), the 
particle density in a pure pair plasma. If we compute 
Lkin by using equation (46) and (50), we obtain the 
following ratio 

AT* (pair) _ L kin (pair) TOpC 2 + (7-(nml))m e c 2 
AT*(nml) Lkin(nml) 2(7_ (pair)) m c 2 

= I^^il, (52) 
7min(pair) 

where 7 m i n (pair) and 7 m i n (nml) refer to 7 m i n in a 
pair and a normal plasma, respectively; (7- (pair)) 
and (7_(nml)) refer to the averaged Lorentz fac- 
tor of electrons in a pair and a normal plasma, re- 
spectively. Therefore, if 7min(pair) is comparable to 
7 m in(nml) <~ 100, there is little difference between 
./V*(pair) and JV*(nml). 

3.6. Summary of the Method 

Let us summarize the main points that have been 
made in §§ 2 and 3. 

(1) We can compute the magnetic field strength (i.e., 
line B in fig. 1) from the surface brightness condition. 

(2) We can also constrain the proper electron density 
and the magnetic field strength (i.e., line C) from the 
synchrotron self-absorption constraint. 

(3) Combing (1) and (2), we obtain the electron den- 
sity, A"* (SSA), for each radio-emitting component. 

(4) Observing the core size, ^d,corc( zy )i at a single fre- 
quency, v, and substituting equation (50) into equa- 
tion (46), we obtain the kinetic luminosity, L^ n . 

(5) Assuming a normal-plasma dominance, we obtain 
the proper electron density, A*(nml) (i.e., line D 2 ) 
from Lkin by equation (51). 

(6) If N* (nml) < N* (SSA) holds, we can rule out the 
possibility of a normal-plasma dominance and that of 
a pair-plasma dominance with 7 m ; n <~ 100 or greater. 
That is, A c *(nml) < A C *(SSA) indicates the domi- 
nance of a pair plasma with 7 min <C 100. 



4. Application to the 3C 345 Jet 

Let us apply the method described in the previ- 
ous sections to the radio-emitting components in the 
3C 345 jet and investigate the composition. This 
quasar (z — 0.595; Hewitt & Burbidge 1993) is one of 
the best studied objects showing structural and spec- 
tral variabilities on parsec scales around the compact 
unresolved core (for a review, see e.g., Zensus 1997). 
At this redshift, 1 mas corresponds to 5.85ft -1 pc, 
and 1 mas yr -1 to /3 app = 30.3/i -1 . 

4.1. Kinetic luminosity 

To deduce Lkin, we first consider 0d,corc^ /^o, where 
k r = 1 is adopted in equation (50). From the re- 
ported core size at 22.2 GHz at 6 epochs by Zensus et 
al. (1995), at 5 epochs by Unwin et al. (1997), and 
at 3 epochs by Ros et al. (2000), we can deduce the 
averaged core size, 0d,corc, as 0.296 mas (fig. 4); here, 
we evaluate the core size with for the latter 

three epochs, where a and b refer to the major and 
minor axes at the FWHM of the elliptical Gaussian 
presented in Ros et al. (2000). From #d,corc = 0.296 
mas, we obtain ^d.corc^ /v Q = 6.18 x 10~ 14 as the 
averaged value over the 14 epochs at 22 GHz. In 
the same manner, we obtain the following averaged 
values at different frequencies: 0d,core = 0.258 mas 
and #d,corc*V^o — 3.64 x 10~ 14 over 3 epochs at 15 
GHz (Ros et al. 2000); 6> d , coro = 0.436 mas and 
#d,coiW^o = 4.38 x 10~ 14 over 7 epochs at 10.7 GHz 
(1 epoch from Unwin et al. 1994; 4 epochs from Zen- 
sus et al. 1995; 2 epochs from Gabuzda et al. 1999); 
#d,core = 0.343 mas and 6& tCaie v/vQ = 2.71 x 10~ 14 
over 7 epochs at 8.4 GHz (1 epoch from Unwin et al. 
1994; 3 epochs from Zensus et al. 1995; 3 epochs 
from Ros et al. 2000); #d,core = 0.423 mas and 
^d,corc^/^o = 1-98 x 10~ 14 over 11 epochs at 5.0 GHz 
(4 epochs from Brown et al. 1994; 1 epoch from Un- 
win et al. 1994; 3 epochs from Zensus et al. 1995; 3 
epochs from Ros et al. 2000). Taking a weighted av- 
erage of ^d,coro^/^o over the 42 epochs at the different 
5 frequencies, we obtain 

= 4.01 x 10- 14 . (53) 

We next consider the spectral index that reflects 
the energy distribution of the power-law electrons 
in the core. From infrared to optical observations 
(J,H,K,L bands from Sitko et al. 1982; J,H,K bands 
from Neugebauer et al. 1982; K band from Impey et 
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al. 1982; IRAS 12, 25, 60, 100 microns from Moshir 
ct al. 1990), we obtain 



a = -1.15 



(54) 



for the core. This value is consistent with that 
adopted by Zensus et al. (1995). 

Using equations (53) and (54), we obtain from 
equation (50) 



1 x^r 



2(5-2a)/(7-2a) 



v vq n sin ip 
Moreover, a = —1.15 gives 

-2aK n 



= 4.97 x 10" 21 /i~ 157 . 

(55) 



7rC(a) 



/mm ' 



-4/(7-2a) 



14^-0.43„, . -0.56 
mm 



= 1.31 x 10- i4 if- u - 4J 7l 
We thus obtain 



(56) 



L kin = 4.68 x 10 45 / i - 1 - 57 ( 7+ ) (^) _1 ' 56 ^ 



0.57 



xr(r - 1) 



i 



sin p \ 1 + z 



2.05' 



-0.43 



ergs s 1 (57) 



for a normal plasma, that is, Cki n = k 2 (7 + )to p c 3 / (27 m i„) 
in equation (46). Note that the factor (7+) cancels 
out when we compute A*(nml) by equation (51) and 
hence does not affect the conclusions on the matter 
content in the present paper. However, to examine 
the absolute kinetic luminosity of a normal-plasma- 
dominated jet, we have to independently know (7+), 
because we can constrain only the electron density by 
SSA theory. 

4.2. Kinematics near the core 

To constrain Lkin further, we must consider the 
kinematics of the core and evaluate T and ip. In the 
case of the 3C 345 jet, components are found to be 'ac- 
celerated' as they propagate away from the core. For 
example, for component C4, /3 app h < 8 held when the 
distance from the core was less than 1 mas (Zensus 
et al. 1995); however, it looked to be 'accelerated' 
from this distance (or epoch 1986.0) and attained 
Pupph ~ 45 at epoch 1992.5 (Lobanov & Zensus 1994). 
It is, however, unlikely that the fluid bulk Lorentz fac- 
tor increases so much during the propagation. If /3 app 
reflected the fluid bulk motion of C4, the bulk Lorentz 
factor, T, had to be at least 45 until 1992.5. Then 



/3 a pp ~ 6 during 1981-1985 indicates that the view- 
ing angle should be less than 0.05°, which is probably 
much less than the opening angle. Analogous 'accel- 
eration' was reported also for C3 (Zensus et al. 1995) 
and C7 (Unwin ct al. 1997). 

On these grounds, we consider the apparent 'ac- 
celeration' of the components in their later stage of 
propagation on VLBI scales are due to geometrical ef- 
fects such as 'scissors effects' (Hardee & Norman 1989; 
Fraix-Burnet 1990), which might be created at the in- 
tersection of a pair of shock waves. Such shock waves 
may be produced, for example, by the interaction be- 
tween relativistic pair-plasma beam and the ambient, 
non-relativistic normal-plasma wind (Sol, Pelletier & 
Asseo 1989; Pelletier & Roland 1989; Pelletier & Sol 
1992; Despringre & Fraix-Burnet 1997; Pelletier & 
Marcowith 1998). In their two-fluid model, the pair 
beam could be destructed at a certain distance from 
the core on VLBI scales due to the generation of Lang- 
muir waves. 

Then, how should we constrain T and ip for the 
3C 345 jet? In this paper, we assume that the radio- 
emitting components before their rapid acceleration 
represent (or mimic) the fluid bulk motion. Under 
this assumption, Stcffcn et al. (1995) fitted the mo- 
tions of C4 (during 1980-1986) and C5 (during 1983- 
1989) when the components are within 2 mas from the 
core by a helical model. They found T = 5.8 for C4 
and 4.6 for C5 with ip = 6.8°, assuming h — 1.54 (i.e., 
.Ho = lOOkm/s/Mpc). These Lorentz factors will be 
<~ 9 and <~ 7 for C4 and C5, respectively, if ft = 1 
(i.e., Hq = 65km/s/Mpc). Subsequently, Qian et al. 
(1996) investigated the intrinsic evolution of C4 under 
the kinematics of T = 5.6 and ip = 2° - 8°, assuming 
h= 1.54; the range of viewing angles is in good agree- 
ment with ip = 2° (for C4) - 4° (for C2) at 1982.0 
obtained by Zensus et al. (1995), who used a larger 
value of T = 10 under h = 1.54, which corresponds to 
r ~ 15 for ft = 1. 

Unless a substantial deceleration takes place, the 
assumption of a constant Lorentz factor is justified 
as the first order of approximation. In the last part 
of § 4.3, we will examine the case when the Lorentz 
factor is different from this assumed value and con- 
sider how the conclusion of the composition depends 
on T. Close to the core, we adopt <p = 2°, because the 
viewing angle is suggested to decrease with decreas- 
ing distance from the core (Zensus et al. 1995; Unwin 
et al. 1997). If p is less than 2° in the core, Lkin, and 
hence N* (nml) further decreases; that is, the normal- 
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plasma dominance can be further ruled out. In short, 
we apply K = 0.1, T = 15 and <p> = 2° in equa- 
tion (57) . Assuming a normal-plasma dominance, we 
adopt 7 m ; n = 10 2 for electron energy distribution. 

4.3. Composition of individual components 

We can now investigate the composition of the 
radio-emitting components in the 3C 345 jet, using 
their spectral information. We assume that the jet 
is neither accelerated nor decelerated and apply the 
Lorentz factor obtained for the (unresolved) core to all 
the (resolved) jet components. To compute A* (SSA), 
we further need ip for each component at each epoch. 
(Note that p = 2° is assumed for the core, not for 
the components.) For this purpose, we assume that 
/?app reflects the fluid motion if /3 app < Vr 2 — 1 and 
that <f is fixed after /3 app exceeds \/r 2 — 1. That is, 
we compute p from 



tantp = < 



2/3 app 

/3app 2 + 6* - 1 



I Vr 2- T 



if/? app < Vr^T 



if /3 app > Vr^i 



(58) 



The constancy of p when /3 app > Vr 2 — 1 = 15.0 
could be justified by the helical model of Stcffen et 
al. (1995), who revealed that the viewing angle does 
not vary significantly after the component propagates 
a certain distance (about 1 mas in the cases of com- 
ponents C4 and C5) from the core. As we will see 
in table 2 (later in this subsection), tp increases with 
increasing distance from the core and do not deviate 
significantly from the assumed value for the core (2 
degree)." 

The synchrotron self-absorption spectra of individ- 
ual components are reported in several papers. They 
were first reported in Unwin et al. (1994), who pre- 
sented (i> m ,S m ,a) and 0d (or £ in their notation) of 
components C4 and C5 at epoch 1990.7. Subse- 
quently, Zensus et al. (1995) gave (v m ,S m ,a), 0d, 
and /3 app /i for C2, C3, and C4 at 1982.0. However, 
the errors are given only for 0d; therefore, we eval- 
uate the errors of the output parameters, S, tp, B, 
A* (nml), and A* (nml) /AT* (SSA), using the errors in 
d alone for C2, C3, and C4 at 1982.0. Later, Unwin 
et al. (1997) presented (i/ m ,S m ,a), 6<\, and j3 apv h of 
C5 (at 1990.55) and C7 (at 1992.05, 1992.67, 1993.19, 
and 1993.55). More recently, Lobanov and Zensus 
(1999) presented a comprehensive data of {v m ,S m ,a) 
of C3, C4, and C5 at various epochs. We utilize 



d = 0.114 + 0.0658(epoch - 1979.50) for C4 (Qian 
et al. 1996), d = 0.114 + 0.0645(epoch - 1980.00) 
for C5, which is obtained from the data given in 
Zensus et al. (1995) and Unwin et al. (1997), and 
d = 0.09 + 0.45(psin<p/mas) for C3 (Biretta et al. 
1986), where psimp refers to the projected distance 
from the core. 

The results for individual components at various 
epochs are given in table 2 for K = 0.1. It follows 
that the ratio, A* (nml) /A* (SSA) becomes less than 
1 if the magnetic field is an acceptable strength, (e.g., 
B < lOOmG). For C4, we cannot rule out the possi- 
bility of a normal-plasma dominance at epochs 1985.8 
and 1988.2. Nevertheless, at these two epochs, un- 
natural input parameters (a = —0.2 at epoch 1985.8 
and S m — 0.8 Jy at epoch 1988.2) give too large 
B (> 380mG) for a moderate <5(< 30). If we in- 
stead interpolated 5 m and a at epochs 1985.8 and 
1988.2 from those at epochs 1983.4 and 1990.7, we 
would obtain reasonable values B = 70 mG and 
A *(nml)/A C *(SSA) = 0.28 at epoch 1985.8, and B = 
50 mG and N* (nml) /A* (SSA) = 0.046 at epoch 
1988.2. On the other hand, C5 at 1984.2 gives un- 
usual value A* (SSA) = 3.7 x 10 5 cnr 3 . It is due to 
the unnaturally small value of i / m (= 2.5 GHz) in the 
early stage of its evolution. 

For an assumed Lorentz factor r — 15 throughout 
the jet, we present the results of N* (nml)/A* (SSA) as 
a function of the projected distance psimp in figure 5. 
The horizontal solid line represents N* — A* (nml). 
It should be noted that (7_) ~ 10 3 holds in equa- 
tion (51) for a normal plasma. Thus, the solid line 
gives the upper limit of N* (nml)/A* (SSA) . What has 
be noticed is that the possibility of a normal-plasma 
dominance can be ruled out if A* (nml) /A* (SSA) <C 1 
holds. There is also depicted a horizontal dashed line 
representing 0.01 A* (nml) /A* (SSA) w A c *(pair)/A C *(SSA) 
(see eq. [52]). If the ratio appears near the dashed 
line, it indicates that the homogeneous component 
is pair dominated. The horizontal dotted line cor- 
responds to (2 / 1836) A* (nml) /N* (SSA) . If the ratio 
appears under this line, it means that Lki n is underes- 
timated for that particular component. If z/ m resides 
within the observed frequency range, the point and 
error bar is indicated by thick pen, whereas if v m re- 
sides outside of the observed frequency range (i.e., if 
v-ca is extrapolated and hence contains a large error 
together with S m ) it is indicated by thin one. The 
two points appearing above the solid line correspond 
to C4 at 1985.8 and 1988.2, which have the unnat- 
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ural input parameters as discussed in the foregoing 
paragraph. 

It follows from the figure that the upper bound 
of the 68%-error bars for the 11 epochs appear be- 
low 0.1, and hence that the jet is likely pair-plasma 
dominated, provided T ~ 15. (C5 at 1984.2 is 
not counted in the 11 epochs, because it gives too 
small N* (nml) /N* (SSA) value.) For a smaller T, 
7V*(nml)/iV*(SSA) further decreases. However, in 
this case, the ratio appears lower than the dotted line, 
indicating that ikin is underestimated for those jet 
components with T < 15, or that (7+) 3> 1 holds for 
a normal plasma composition. The ratio shows no 
evidence of evolution as a function of the projected 
distance, p sin ip, under the assumption of a constant 

r. 

It is noteworthy that N* (nml) /N* (SSA) is pro- 
portional to (imin/lOO) -2 " -1 ' 56 for the spectral in- 
dex of —1.15 for the core. Therefore, for a typical 
spectral index a ~ —0.75 for the jet components, the 
dependence on 7 m ; n virtually vanishes. That is, the 
conclusion does not depend on the assumed value of 
7min ~ 100 for a normal plasma. 

If we assume instead T = 20 (rather than T = 15) 
and p — 2° , we obtain the following kinetic luminosity 



L kin = 1.3xl0 46 /i 



-1.57 



/ x / 7min \ 

(7+> lloo ) 



-1.56 



K 



0.56, 



ergs s 
(59) 

Evaluating the viewing angles of individual compo- 
nents by equation (58) with T = 20, we can compute 
N* (nml) /N* (SSA) as presented in figure 6. It follows 
from this figure that we cannot rule out the possibility 
of a normal-plasma dominance for such large Lorentz 
factors. The greater T we assume, the greater be- 
comes A c *(nml)/A C *(SSA). 

5. Application to the 3C 279 Jet 

Let us next consider the 3C 279 jet. This quasar 
(z = 0.538) was the first source found to show su- 
perluminal motion (Cotton ct al. 1979; Unwin et al. 
1989). At this redshift, 1 mas yr _1 corresponds to 
/3 app = 28.2/j- 1 . 

5.1. Kinetic luminosity 

In the same manner as in 3C 345 case, we first 
consider #d,corc*V^o- Since the turnover frequency 
of the core is reported to be about 13 GHz (Un- 
win ct al. 1989), we measure #d,core at v < 10.7 



Core size of 3C345 jet at 22 GHz 
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Fig. 4. — Angular core size of the 3C 345 core at 22.2 
GHz. The horizontal line represents the bet fit value. 
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distance from core [mas] 

Fig. 5.— The ratio N* (nml) /N* (SSA) as a function 
of the distance from the core. Above the solid, hori- 
zontal line, the dominance of a normal plasma is al- 
lowed. T = 15 is assumed. 
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GHz (i.e., optically thick frequency range). At 10.7 
GHz, #d,corc is reported to be 0.79 mas from five- 
epoch observations by Unwin et al. (1989), and 0.61 
mas from three-epoch observations by Carrara et al. 
(1993). With one-epoch VSOP (VLBI Space Obser- 
vatory Programme) observation, #d,corc = 1-8 x 0.28 
mas and #d,corc = 1-8 x 1.95 mas were obtained at 4.8 
GHz and 1.6 GHz, respectively (Pincr et al. 2000). 
As the weighted mean over the 10 epochs, we obtain 



6.6 x 10" 



(60) 



which leads to 
1 x^r 



v n sin ip 



3.7 x 10 



^d,corc 

mas 

13 



V 

v 



(61) 



We next consider the spectral index of the core in 
the optical thin frequencies. From mm to IRAS (90- 
150 GHz and 270-370 GHz) observations, Grandi et 
al. (1996) reported 



a = -1.2. 



(62) 



Adopting this value as the spectral index that reflects 
the energy distribution of the power-law electrons in 
the core, we obtain 



1 X^rv 

i>q r\ sin ip 



2(5-2a)/(7-2ce) 



= 2.7 x 10 _20 /i 



-1.57 



(63) 



and 



nC(a) 



-2aK n 

7min 2a+1 n> 



-i -4/(7-2a) 



1.81 x 10~ 14 K 



-0.43 -0.56 
/min 



(64) 



If we assume a normal plasma dominance, we obtain 

-1.56 



V ioo / 

x r ( r - !) 



57 



sin ip \ 1 + z 



5.2. Kinematics near the core 



0.43 



ergs s . 
(65) 



We now consider T and if near to the core to con- 
strain Lkin further. The apparent superluminal veloc- 
ities of the components in the 3C 279 jet do not show 



the evidence of acceleration and are roughly constant 
during the propagation. For example, Carrara et al. 
(1993) presented /? app = 4.51/j" 1 for C3 during 1981 
and 1990 over 15 epochs at 5, 11, and 22 GHz and 
/3 app = 4.23/j- 1 for C4 during 1985 and 1990 over 
15 epochs at 11 and 22 GHz. We thus assume that 
the apparent motion of the components represent the 
fluid bulk motion and that both T and if are kept 
constant throughout this straight jet. Since /3 app does 
not differ very much between C3 and C4, we adopt 
/3 a pp = 4.51 as the common value and apply the same 
T and ip for these two components. The errors in- 
curred by the difference in /3 app are small. 

The lower bound of T can be obtained as T > 
\J Aipp 2 + 1 = 4.61. The upper bound of T, on the 
other hand, can be obtained by the upper bound of 6 
for a given (3 app = 4.51. In general, the upper bound 
of S is difficult to infer. In the case of 3C 279, however, 
it is reasonable to suppose that a substantial fraction 
of the X-ray flux in the flaring state is produced via 
synchrotron self Compton (SSC) process. Therefore, 
the Doppler factors <~ 3.9 (Mattox et al. 1993) and 
- 5 (Henri et al. 1993) derived for the 1991 flare, 
and ~ 6.3 (Wehrle et al. 1998) for the X-ray flare in 
1996, give good estimates. We thus consider that 8 
docs not greatly exceed 10 for the 3C 279 jet. Im- 
posing S < 10, we obtain T < 6.06 from /3 app = 4.61. 
Combining with the lower bound, we can constrain 
the Lorentz factor in the range 4.6 < T < 6.1. 

As case 1, we adopt the smallest bulk Lorentz 
factor, T = 4.61, which results in ip = 12.5° and 
5 = 4.6 for /3 app = 4.51. It is worth comparing this 
5 with that measured in the quiescent state; Unwin 
derived 5 = 3.6^'° so that the calculated X-ray flux 
may agree with the observation. Moreover, from the 
RXTE observation in the 0.1-2.0 keV band during 
1996.06 and 1996.11, Lawson and M c Hardy (1998) 
obtained 2/iJy in the quiescent state during 1996.49- 
1996.54 from 16-epoch observations. Applying this 
X-ray flux density to the spectral information (v m , 
S m , a) and #d (table 3) for component C4 at epochs 
1987.40 and 1989.26 and component C3 at 1984.05, 
we obtain S > 1.9, > 4.0, and > 3.3, respectively. 
Therefore, we can regard that the value of 5 = 4.6 are 
consistent with the X-ray observations. In this case, 
equation (65) gives 

i ki n = 7.2xl0 46 ^- 57 ( 7+ )(^) _1 - 6 

(66) 

As case 2, we adopt a large Lorentz factor of 



h 16 ergs s 1 (case 1) 
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r = 6.0. In this case, /3 a pp — 4.51 gives 5 = 9.8 and 
ip = 4.45°. There is, in fact, another branch of solu- 
tion that gives smaller Doppler factor of 2.2; however, 
this solution gives too large AT* (SSA) 's (~ 10 5 cm~ 3 ). 
We thus consider only the solution giving 5 = 9.8 as 
case 2, which results in 

L kin = 3.5xl0 46 A°- 57 < 7+ ) (^"'V^ergs s"* (, 

(67) 

We will examine these two cases in the next subsec- 
tion. 

5.3. Composition of individual components 

For the jet component C3, Unwin et al. (1989) 
presented v m = 6.8 GHz, S m = 9.4 Jy, a = —1.0, 
and #d = 0.95 mas at epoch 1983.1 (table 3). It was 
also reported in their paper that the flux densities 
were 5.13 ± 0.14 Jy at 5.0 GHz (at epoch 1984.25), 
4.49 ±0.17 Jy at 10.7 GHz (1984.10), and 2.58 ±0.18 
Jy at 22.2 GHz (1984.09) for C3. We model-fit the 
three-frequency data by the function 

S v = A l [1 - cxp(-A 2 ^- 2 - 5 )] . (68) 

Assuming a = —1.0, wc obtain ^ m = 6.6 GHz and 
S m = 5.9 Jy as the best fit at epoch 1984.10. 

For C4, Carrara et al. (1993) presented is m ~ 11 
GHz, S m <~ 4.3 Jy, a = —0.9, and 8d ~ 0.6 mas from 
their 1989-1990 maps at 5, 11, and 22 GHz. At epoch 
1987.4, we use the flux density of 1.43 ± 0.17 Jy at 22 
GHz (Carrara et al. 1993) and that of 3.95 ± 0.20 Jy 
at 5 GHz (Gabuzda et al. 1999). We extrapolate the 
flux density at 11 GHz from those at 1988.17, 1989.26, 
and 1990.17 presented in Carrara et al. (1993) to ob- 
tain 3.60±0.20 Jy. Assuming a — —0.9, we can model 
fit the three-frequency data to obtain v m = 6.4 GHz 
and S m = 4.4 Jy. At this epoch (1987.4), component 
C4 is located about 1 mas from the core; therefore, 
the angular size vs. distance relation (Carrara et al. 
1993) gives 6d ~ 0.6 mas. The input parameters for 
C3 and C4 at these 4 epochs are summarized in ta- 
ble 3. 

In case 1 (see § 5.2), we obtain AT* (SSA) = 4.3 x 
10 3 and 3.0 x 10 4 , for C3 at 1983.1 and 1984.10, and 
8.0 x 10 3 and 6.4 x 10 1 for C4 at 1987.4 and 1989.26, 
respectively. These large values of A* (SSA) in the 
first three epochs result in such (un-physically) small 
values of AT* (nml)/iV* (SSA) as 1.4 x 10~ 4 , 2.9 x 10~ 5 , 
and 1.9 x 10 -4 . Since even a pair-plasma dominance 
should be ruled out for N* (nml) /AT* (SSA) < 100 (see 



eq. [52]), we consider that the value of 6(= 4.6) is 
underestimated. 

We next examine case 2, in which a larger Doppler 
factor (= 9.8) is adopted. In this case, we obtain 
reasonable values of B and A* (SSA) as presented in 
table 3. It follows that the ratio A c *(nml)/A C *(SSA) 
becomes 4.3 x 10~ 3 and 8.7 x 10~ 3 , for C3 at 1983.1 
2). and 1984.10, and 5.0 x 10~ 3 and 0.62 for C4 at 1987.4 
and 1989.26, respectively. If 5 exceeds 10 (or equiv- 
alently, if T exceeds 6), the density ratio for C4 at 
1989.26 becomes close to unity; therefore, we cannot 
rule out the possibility of a normal-plasma dominance 
for such large Doppler factors. 

On these grounds, we can conclude that the jet 
components are dominated by a pair plasma with 
7min -C 100, provided that 5 < 10 holds, in the first 
three epochs. For the last epoch (1989.26), the 68%- 
error bar of N* (nml) /N* (SSA) appears above 0.1; 
thus, we cannot rule out the possibility of a normal- 
plasma dominance. If 5 is as small as 5, not only 
a normal-plasma dominance but also a pair-plasma 
dominance should be ruled out for the first three 
epochs. Thus, we consider S is greater than 5 for 
the 3C 379 jet. 

6. Discussion 

In summary, we derived a general scheme to infer 
the kinetic luminosity of an AGN jet using a core 
size observed at a single VLBI frequency. The kinetic 
luminosity gives an electron density as a function of 
the composition of the jet. If the density deduced 
under the assumption of a normal-plasma dominance 
becomes much less than that obtained independently 
from the theory of synchrotron self-absorption, we can 
exclude the possibility of a normal-plasma dominance. 
Applying this method to the 3C 345 jet, we found that 
components C2, C3, C4, C5, and C7 are dominated 
by a pair plasma with 7 m ; n <C 100 at 11 epochs out of 
the total 21 epochs examined, provided that T < 15 
holds throughout the jet. We also investigated the 
3C 279 jet and found that components C3 and C4 are 
dominated by a pair plasma with 7 m i n <C 100 at three 
epochs out of the four epochs examined, provided that 
5 < 10. 

It is noteworthy that the kinetic luminosity com- 
puted from equations (46) and (50) has, in fact, weak 
dependence on the composition. Substituting Cki n for 
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a pair and a normal plasma, we obtain 

L kin (pair) m c (t_ (pair)) / 7 min (pair) 

^kin(nml) m p /27 min (nml) 

1 (7- (pair)) 7 min (nml) 
9.18 7min (pair) 100 



.(69) 



It follows from equation (40) that (7. (pair) )/7 m in (pair) 
3 holds for a — —0.75, for instance. Thus, we ob- 
tain comparable kinetic luminosities irrespectively of 
the composition assumed, if we evaluate them by the 
method described in § 3. 

Let us examine the assumption of K ~ 0.1. For a 
typical value of a = —0.75, we obtain ((7_)/7 m i n ).fir = 
(<5 cq /<5) 17 / 2 , where 8 cq refers to the "equipartition 
Doppler factor" defined by (Readhead 1994) 

S cq = |10- 3 - 6 " F(a) 3i d(z) 2 



x (1 + z) 



17-2a,, -35-2a c 16/) -34 "1 l/(13-2a) 



where 



d(z) 



0.65/i 



qoz + (q - 1)(-1 + ^2q z + 1) ' 



(70) 



(71) 



and F(a) is given in Scott and Readhead (1977); f m , 
S m , and 6<\ are measured in GHz, Jy, and mas, re- 
spectively. Therefore, we can expect a rough energy 
equipartition between the radiating particles and the 
magnetic field if <5 oq is close the 5 derived indepen- 
dently. 

Substituting the input parameters presented in ta- 
ble 2, we can compute 8 cq s for the 3C 345 jet. The 
computed <5 cq ranges between 1.1 and 67 except for 
component C5 at 1984.2, which gives <5 oq = 736 due 
to unnaturally small v m . If we singled out this epoch, 
we would obtain <5 eq = 23 ±9. The resultant magnetic 
field strength, B = 6.7 ± 6.4 mG, is reasonable. For 
the 3C 279 jet, we obtain S cq = 9.5 ± 1.5 and 12 ± 8 
mG for the four epochs examined. 

Since <5 cq = 23 ± 9 for 3C 345 is consistent with 
those presented in table 2 and since <5 oq = 9.5 ± 1.5 for 
3C 279 is consistent with those in table 3, we can ex- 
pect a rough energy equipartition in these two blazer 
jets. However, it may be worth noting that the en- 
ergy density of the radiating particles dominates that 
of the magnetic field if 5 becomes much smaller than 
<5 oq . Such a case (5 = 4 ~ 12) was discussed by Unwin 
et al. (1992; 1994; 1997) for the 3C 345 jet, assuming 
that a significant fraction of the X-rays was emitted 



via SSC from the radio-emitting components consid- 
ered. 

We finally discuss the composition variation along 
the jet. In the application to the two blazers, we as- 
sumed constant Lorentz factors (r) throughout the 
jet. However, T may in general decrease with increas- 
ing distance from the central engine. If the decrease 
of T is caused by an entrainment of the ambient mat- 
ter (e.g., disk wind) consisting of a normal plasma, it 
follows from equation (51) that A^*(nml) and hence 
TV* (nml)/N* (SSA) increases at the place where the 
entrainment occurs. (Note that the kinetic luminosity 
will not decrease due to an entrainment in a station- 
ary jet.) Therefore, the change of composition from a 
pair plasma into a normal plasma along the jet may 
be found by the method described in this paper. 

To study further the change of the composition, we 
need detailed spectral information (f m , 5 m , a) of in- 
dividual components along the jet. To constrain the 
spectral turnover accurately, we must model-fit and 
decompose the VLBI images from high to low fre- 
quencies; in another word, spatial resolution at low 
frequencies is crucial. Therefore, simultaneous obser- 
vations with ground telescopes at higher frequencies 
and with space+ground telescopes at lower frequen- 
cies are essential for this study. 
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Fig. 6. — The same figure as figure 3, but the different 
Lorentz factor T = 20 is adopted. 
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Table 2: Input parameters and the results for the 3C 345 jet when T = 15. 



epoch 


p sin ip' 




"1 




a* 




6 




B 


jV e *(SSA) 


K(nml)/AT e *(SSA) 




mas 




GHz 


Jy 




mas 




deg 


mG 


cm 3 


(68 % errors) 


-xf 












component 


C2 










1982.00 a 


4.9 


12.9 


1.5 


2.0 


-0.60 


2.15 ± .18 


22.5 


2.2 


1.4E1 


6.2E-4 


2.5E-1 ( 1.3E-1, 4.5E-1) 














component 


C3 










^wltso 15 


2.0 


8.7 


2.7 ±0.3 


2.4 ±0.2 


-0.80 


1.00 


27.2 


1.2 


9.7E0 


2.9E-2 


2.4E-2 ( 8.6E-3, 6.1E-2) 


1^82.00 a 


2.2 


9.2 


2.6 


2.1 


-0.70 


0.97 ± .07 


26.8 


1.3 


9.5E0 


1.5E-2 


5.1E-2 ( 2.8E-2, 8.7E-2) 














component 


C4 










4©82.00 a 


0.4 


6.2 


14.6 


7.6 


-0.30 


0.29 ± .02 


28.6 


.82 


2.7E1 


2.2E-2 


3.9E-1 ( 2.5E-1, 5.8E-1) 


<^£83.40 b 


0.7 


6.0 


11.6 ±0.4 


7.5 ±0.6 


-0.50 


0.37 


28.7 


.81 


2.9E1 


4.0E-2 


1.3E-1 ( 8.6E-2, 1.9E-1) 


1985.80 b 


1.0 


7.9 


9.4 ±0.6 


1.9 ±0.1 


-0.20 


0.53 


27.7 


1.1 


3.9E2 


2.3E-5 


1.1E2 ( 7.2E1, 1.6E2) 


r+968.20 a 


2.0 


14.2 


6.5 ± 1.4 


0.8 ±0.1 


-0.60 


0.69 


19.7 


2.8 


1.2E3 


2.4E-5 


6.1E1 ( 9.5E0, 2.8E2) 


t&90.70 c 


4.3 


14.0 


2.0 ±0.5 


2.0 ±0.5 


-0.75 ±.15 


0.85 


20.3 


2.6 


1.2E0 


3.8E0 


2.6E-4 ( 1.4E-5, 3.8E-3) 


tl992.50 b 


4.3 


13.4 


2.7 ± 1.6 


1.6 ±0.1 


-0.70 


0.97 


21.7 


2.4 


1.6E1 


1.5E-2 


5.0E-2 ( 1.5E-3, 5.5E0) 



component C5 



T^4.20 b 


0.3 


5.0 


2.5 ± 0.4 


4.9 ± 1.1 


-0.80 


0.27 


29.1 


.67 


9.0E-3 


3.7E+5 


2.6E-8 


; 4.8E-9, 1.4E-7) 


1985.80 b 


0.6 


6.4 


6.0 ±0.3 


3.3 ±0.3 


-0.60 


0.37 


28.5 


.86 


5.7E0 


4.0E-1 


1.3E-2 


; 7.6E-3, 2.2E-2) 


4J87.30 b 


0.9 


7.5 


7.0 ± 0.6 


5.1 ±0.2 


-0.80 


0.47 


28.0 


1.0 


1.3E1 


6.2E-1 


5.2E-3 


; 2.4E-3, 1.0E-2) 


CT^8.20 b 


1.1 


8.0 


7.5 ±0.4 


5.6 ±0.5 


-0.70 


0.53 


27.6 


1.1 


2.4E1 


7.9E-2 


3.2E-2 


; 1.8E-2, 5.6E-2) 


^2£9.20 b 


1.4 


8.6 


2.5 ± 1.4 


4.7 ±0.6 


-0.90 


0.59 


27.2 


1.2 


2.0E-1 


9.1E2 


2.3E-6 


; 6.1E-8, 3.5E-4) 


T5&0.30 b 


1.7 


9.3 


2.5 ± 1.2 


3.6 ±0.1 


-0.90 


0.66 


26.7 


1.3 


5.2E-1 


9.3E1 


1.8E-5 


; 5.0E-7, 1.1E-3) 


ie9u.55 d 


1.8 


9.5 


2.7 ±0.5 


3.2 ±0.5 


-0.75 


0.80 


26.5 


1.4 


2.2E0 


8.2E-1 


1.4E-3 


; 2.4E-4, 6.5E-3) 


Seo.70 c 


1.7 


11.3 


2.7 ± 1.5 


3.2 ±0.9 


-0.75 ±.15 


0.80 ± .20 


24.8 


1.7 


2.1E0 


1.4E0 


7.9E-4 


; 1.5E-5, 4.2E-1) 


iS92.50 b 


2.3 


11.2 


2.2 ± 0.8 


1.5 ±0.3 


-0.90 


0.81 


24.9 


1.7 


3.4E0 


8.5E-1 


1.3E-3 


; 4.7E-5, 3.3E-2) 


C/J . 












component 


C7 














0.1 


2.8 


12.8 ±0.5 


4.6 ±0.5 


-0.75 


0.20 ± .04 


29.7 


.36 


1.1E1 


2.7E0 


6.6E-3 


; l.OE-3, 3.1E-2) 




0.2 


6.0 


12.5 ± 1.0 


7.0 ±0.5 


-0.75 


0.35 ± .02 


28.7 


.80 


3.9E1 


2.1E-1 


2.8E-2 


; 1.1E-2, 6.1E-2) 


|jae3.19 d 


0.4 


10.5 


11.6 ±0.5 


5.1 ±0.5 


-0.75 


0.41 ± .02 


25.7 


1.6 


8.6E1 


6.0E-2 


7.0E-2 


; 3.5E-2, 1.3E-1) 


^3.55 d 


0.5 


14.5 


11.0 ± 1.5 


3.1 ± 1.0 


-0.75 


0.38 ± .02 


20.9 


2.6 


1.1E2 


8.8E-2 


5.3E-2 


; 8.9E-3, 4.1E-1) 



^or the input parameters, 1 a errors are presented if they are given in the litearture. 

a Input parameters are cited from Zensus et al. (1995). 

b Input parameters are cited from Lobanov and Zensus (1999). 

c Input parameters are cited from Unwin et al. (1994). 

d Input parameters are cited from Unwin et al. (1997). 
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Tabic 3: Input parameters and the results for the 3C 279 jet in case 2. 



epoch 


psiiiip 






5 ni 


a 


o d 


6 


f 


B 


7V C *(SSA) 


7V *(nml)/7V c *(SSA) 




mas 




GHz 


Jy 




mas 




deg 


mG 


cm~ 3 


(68 % errors) 














component 


C3 










1983.10 a 


1.4 


4.5 


6.8 


9.4 


-1.0 


0.95 


9.8 


4.5 


1.8E1 


3.5E1 


4.3E-3 ( 3.6E-3, 5.0E-3) 


1984.10 b 


1.4 


4.5 


6.5 


8.9 


-1.0 


0.80 


9.8 


4.5 


8.2E0 


2.5E2 


8.7E-4 ( 7.3E-4, 1.0E-3) 














component 


C4 










1987.40 b 


1.0 


4.2 


6.4 


4.4 


-0.9 


0.60 


9.8 


4.5 


9.8E0 


7.7E1 


5.0E-3 ( 4.1E-3, 5.7E-3) 


1989. 26 c 


1.2 


4.2 


11.0 


4.3 


-0.9 


0.60 


9.8 


4.5 


1.6E2 


6.2E-1 


6.2E-1 ( 2.9E-1, 1.3E0) 



b 



Input parameters are cited from Unwin et al. (1989). 
Input parameters are computed in this paper (see text). 
Input parameters are cited from Carrara et al. (1993). 
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